Analysis of three-neutrino oscillations in the full mixing angle space 
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We use a model, with no CP violation, of the world's neutrino oscillation data, excluding the 
LSND experiments, and search the full parameter space (0 < 6*12 < vr/2; — n/2 < 813 < 7r/2; and 
< 623 < 7r/2) for the best fit values of the mixing angles and mass-squared differences. We find 
that the mixing angle #13 is bounded by —0.15 < #13 < 0.20 with an absolute minimum at 813 — 0.12 
and a local minimum at —0.04. The importance of the negative 813 region and this structure in the 
chi-square space has heretofore been overlooked because the factorization approximation commonly 
employed yields oscillation probabilities that are a function of sin 2 813. 
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The observation of neutrino oscillations requires a fun- 
damental modification of the electroweak theory. The 
simplest, but not totally consistent, method for accom- 
modating neutrino oscillations into the theory is to in- 



troduce a posteriori a mass matrix and unitary mixing 
matrix. The standard representation of the three neu- 
trino mixing matrix is 



C12C13 S12C13 si3e~ lS 

Uak — > ( -S12C23 - Ci 2 S23Sl3e lS Cl 2 C 2 3 - Sl 2 S23Sl3e 1 ' 5 S23C13 ) , (1) 
S12S23 - Ci2C23Sl 3 e Jl5 -C12S23 - Si 2 C23Si 3 e iS C23C13 



where c a k — cos6 a k, s a k — $mQ a k, and 8 is the CP 
violating phase with 9 a k and 5 real. We order the mass 
eigenstates by increasing mass, and the flavor eigenstates 
are ordered electron, mu, tau. The bounds on the mixing 
angles 6jk are < Qjk < 7r/2 and < 6 < 2tt. In the 
absence of CP violation, the range of the mixing angles 
2] is < Ojk < 7r/2 with 5 — and S — 7r; or equivalently 
3 take only 6 = with < 6 12 < tt/2, -tt/2 < 13 < 
7r/2, and < 6*23 < 7r/2. Experiments find that O13 is 
near zero. The second option produces one contiguous 
allowed region in the parameter space; the former gives 
two disconnected regions for the allowed parameters. In 
particular, oscillation probabilities for 613 negative are 
not related to those for #13 positive. Parameterization of 
oscillation solutions by sin 2 $13 is thus inadequate. 

In vacuo, the probability that a neutrino with energy 
E and flavor a will be detected a distance L away as a 
neutrino of flavor (3 is given by 



V a p{L/E) = 5 aP 



UaiUffjUakUffkWl 2 ^, (2) 



= 1.27Am'j k L/E with L/E expressed in 
1? — mi in units of eV 2 . 



in which </>°| c 

units of m/MeV and Am 2 fc = — uL k 
This probability is then to be integrated over the energy 



spectrum of the neutrinos for each experiment. 

We construct a model of the data and then, within the 
model, look for best fit oscillation parameters throughout 
the full range of permitted mixing angles. Included in 
the model are data for neutrinos from the sun E L I M EIB, 
reactor neutrinos [H |p , atmospheric neutrinos |lfj|. and 
beam-stop neutrino s [III . We, like others, omit from the 
analysis the LSND [12| and Karmenllal experiments. 

Experiments for solar neutrinos |j Li Q historically 
measured the survival probability of electron neutrinos, 
V ee ■ Recent experiments 0, measure two different neu- 
trino interactions which then allow the extraction of "P e e 
and the total solar neutrino flux. The measured total 
is in agreement with the theoretical predictions of the 
standard solar model [14j, We here use the standard 
solar model for the production of neutrinos in the sun. 
Each detector measuring solar neutrinos has a different 
acceptance and thus measures different energy neutrinos. 
In order to reproduce the energy dependence of the sur- 
vival rate of electron neutrinos arriving at the Earth as 
seen in the experiments, we invoke the MSW effect |l5). 
The MSW effect arises because the neutrinos created in 
the sun propagate through a medium with a significant 
electron density. The forward coherent elastic neutrino- 
electron scattering produces an effective change, relative 
to the mu and tau neutrino, in the mass of the electron 
neutrino given by A(r) = \/2 G E p{r) / m n , with p(r) the 



2 



electron density at a radius r, G the weak coupling con- 
stant, and m n the nucleon mass. In the flavor basis, the 
Hamiltonian then becomes 

H mat = UMU* + A(r), (3) 

with M. the (diagonal) mass-squared matrix in the mass 
eigenstate basis and A the 3x3 matrix with the interac- 
tion A(r) as the electron-electron matrix element and ze- 
roes elsewhere. By diagonalizing this Hamiltonian with a 
unitary transformation D a k (r, E) , we define local masses 
and eigenstates as a function or r and E. Care must be 
taken so that D a k(r, E) becomes U a k in the limit of zero 
electron density. In the adiabatic limit, which we use, 
the electron survival probability is 

P e a e d (r 7 E)=J2D(r,E)l k U! k . (4) 

k=l 

Neutrinos are produced throughout the sun by vari- 
ous reactions, each with its own energy spectrum. The 
surviving neutrinos are then detected by detectors which 
have a different acceptance for each energy of the neu- 
trino. The survival probability for an electron neutrino 
in a particular experiment is given by 

Kt = E PT / fi M * / 9i (E) P e a e d (r, E) dE . 

j—l •'O J Ethresh 

(5) 

Here, j labels a particular nuclear reaction; we include 
three reactions - pp, 7 Be, and 8 B. The quantity pj x is the 
probability that in a particular experiment the neutrino 
arose from nuclear reaction j. We take these from the 
analysis of Ref. 01 f° r the solar experiments: chlorine 
0, gallium (Sage,Gallex, and GNO) 0, SNO 0, and 
SNO-salt 0. The function fj(r) is the probability that 
a neutrino is created by reaction j at a radius r Q of the 
sun and is integrated from the center of the sun to the 
solar radius. The function gj(E) is the energy distribu- 
tion of the neutrinos emitted in reaction j. For 7 Be this 
is a delta function at 0.88 MeV; the lower emission line 
does not contribute significantly. For the pp neutrinos, 
the energy distribution times the detector acceptance is 
a relatively narrow function of energy; we set E to its 
average. For 8 B neutrinos, we use the energy distribu- 
tion from the standard solar model |l4| and numerically 
perform the integration. 

For three neutrino mixing, the energy dependence of 
the solar data is well reproduced by the MSW effect with- 
out level crossing. This is true of all the parameter sets 
examined here. The adiabatic approximation is thus jus- 
tified after the fact. 

The reactor experiments that we include are CHOOZ 
and KamLAND ||. KamLAND is unique among re- 
actor experiments as it measures P ee where its predeces- 
sors set limits on I — V ee . It also provides the energy 



Experiment 


Measured 


Ljhj (m/MeVj 


Data 


Chlorine 


Pee 


4.0 x 10 IU 


.337 ± .065 


Gallium 


Pee 


35. x 10 IU 


.550 ± .048 


SNO 


Pee 


2.2 x 10 1U 


.348 ± .073 


SNO-salt 


Pee 


2.2 x 10 IU 


.306 ± .035 


CHOOZ 


Pee 


300. 


> 0.96 


KamLAND 


Pee 


4.1 x 10 4 


.686 ± .006 


K2K 


PflU 


208. 


.55 ±.19 



TABLE I: The experiment, quantity measured, the average 
value of L/E, and experimental data for those quantities fit 
by the model are presented. We have combined the SuperK 
and SNO results into one data point listed as SNO. We list 
Vee for KamLAND but we in fact fit the measured energy 
spectrum. 

spectrum of the neutrinos which tightly constrains the 
small mass squared difference. We list the value of P ee 
in the table, but we actually fit with our model the en- 
ergy spectrum. In order to incorporate the systematic 
error for KamLAND, we introduce a normalization of 
their data, Nxam, and float it constrained by an error 
of six percent. We also include the K2K experiment 
that measures the survival of muon neutrinos over 
a long baseline (250 km) from KEK to the Super-K de- 
tector. The experiment, quantity measured, the value of 
that quantity to which we fit, and the average value of 
L/E for each experiment are given in Table ID 

The Super-Kamiokande experiment |10| has measured 
neutrinos that originate from cosmic rays hitting the 
upper atmosphere. The detector distinguishes between 
e-like (electron and anti-electron) neutrinos and /i-like 
(muon and anti-muon) neutrinos. The rate of e-like neu- 
trinos of energy E arriving at the detector from a source 
a distance L away is 

K e (L, E) = V ee {L, E) + n[E)V^{L, E) , (6) 

and for /i-like neutrinos 

n^L,E)=V^(L,E) + -±-V eii (L,E), (7) 

n(E) 

where n(E) is the ratio of /i-like neutrinos to e-like neu- 
trinos at the source. We incorporate the Super-K atmo- 
spheric neutrinos by utilizing the L / _B-dcpendcncc of lZ e 
and 7?. M given in [Trjj and pictured in Fig. ^ We take 
n(E) to be energy independent and equal to 2.15. As 
the absolute flux of cosmic rays striking the atmosphere 
is not known to within 15%, we introduce as a fit param- 
eter an energy-independent renormalization factor N atm 
that multiplies the experimental values of P, e and 7?. M . 

The ratios P. e and P.^ are convenient for the theorist 
as these are easily calculable. A distinct advantage of 
the atmospheric data is that for the neutrinos arriving 
from directly overhead to those arriving from the oppo- 
site side of the Earth, the value of L/E changes by almost 
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FIG. 1: Data and fits for the Super-K atmospheric exper- 
iment: (a) electron-like and (b) muon-like events. Data is 
represented by points. Fit 1 and fit 2, which is not distin- 
guishable, yield the solid line; solution SS, the dashed-dotted 
line, is only barely distinguishable. 



parameter\fit 


SS 


#1 


#2 


Xdo f 




1.19 


1.21 


012 


0.54 ±0.05 


0.48 


0.49 


013 


< 0.13 


0.12 


-.04 


023 


0.79 ±0.06 


0.81 


0.71 


Am^x x KT^eVf 


8.3 ±0.4 


7.6 


7.7 


Ami 2 x KT^eV) 2 


2.4 ±0.3 


2.6 


2.6 


Natm 




1.00 


1.00 


Nxam 




1.00 


1.00 



TABLE II: The value of Xdof an d the parameters for the 
standard solution (SS) and for the two local minima found 
here. 



four orders of magnitude. This is the only data which 
varies L/E. On the other hand, the source of neutrinos 
from cosmic rays hitting the atmosphere must be mod- 
eled. Also, the relationship between the direction of the 
recoil electrons in the detector and the direction of the 
neutrino initiating the reaction requires additional mod- 
eling. Thus the connection between the quantity mea- 
sured and a simple physical parameterization is indirect 
and difficult to incorporate. The details of the model can 
be found in 

We fit the mixing angles, the mass squared differences, 
N a tm and NKam to the quantities in Tableland to the 
L /E dependence of lZ e and 7\L M pictured in Fig. ^ by 
minimizing chi-squared per degree of freedom, Xdof- ^ n 
Fig. we present A% 2 =: Xdof ~~ Xmin as a function of 
013, where for each value of #13 we have minimized with 
respect to the other parameters. The results are not sym- 




FIG. 2: The value of A\ 2 for a fixed value of 0i3 with all other 
parameters varied. The horizontal lines indicate the 90% and 
95% CL. 



Experiment 


Data 


SS 


#1 


#2 


Chlorine 


.337 ± .065 


.451 


.448 


.454 


Gallium 


.550 ± .048 


.578 


.615 


.623 


SNO 


.348 ± .075 


.395 


.371 


.378 


SNO-salt 


.306 ± .035 


.395 


.371 


.378 


CHOOZ 


> .96 


.98 


.96 


.99 


KamLAND 


.686 ± .006 


.577 


.661 


.670 


K2K 


.55 ±.19 


.60 


.56 


.57 



TABLE III: The experimental results and the predictions for 
each from the models whose parameters are given in Table ITT1 



metric about #13 = 0, and we find two minima. The ab- 
solute minimum is at #13 = 0.12, and there is a second 
local minimum at #13 = —0.04. This asymmetry and 
the existence of the second local minimum would not ex- 
ist if the commonly used factorization approximation to 
the oscillation probabilities were employed, as this gives 
oscillation probabilities that are functions of sin 2 $13. 

In Table |H] we present the oscillation parameters for 
the two minima which we have found and also for a solu- 
tion where Am 2 ^, 612, and #13 are taken from the analysis 
of Ref . 01 j and Am 2 2 and #23 are taken from Ref . 0| . 
We term this latter solution the "standard solution" (SS). 
We see that the parameters we find, particularly for the 
absolute minimum with #13 > 0, are reasonably consis- 
tent with those for the standard solution. We remind the 
reader that we built the model 0] not to extract pre- 
cise values of the oscillation parameters, but to examine 
features of neutrino oscillation phenomenology in a semi- 
quantitative way. Here, we use the model to investigate 
the role of the negative #13 region. 

In Table IIIII we compare the data with the results of 
the fits, and in Fig. we depict the L/E dependence of 
each fit as compared to the atmospheric data. In Fig. 
curves are drawn for all three solutions. However, the 
results are sufficiently similar that the individual curves 
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FIG. 3: Zenith-angle distributions for atmospheric Super-K 
experiment. Points represent the data, and the solid line 
shows the results of Fit 1. (a) Sub-GeV electron- like, (b) 
Sub-GeV muon-like. (c) Multi-GeV electron-like, (d) Multi- 
GeV muon-like and PC events. 



cannot be distinguished. The model treatment of the 
atmospheric data is thus seen to be quite comparable to 
a full analysis. The resulting fits to the data in Table ITTT1 
are also seen to be reasonable. 

In order to further demonstrate that our results are 
reasonable, we calculate the zenith-angle dependence of 
the atmospheric data. Using the ener gies defined for the 
various classes of neutrino events in |10j . we determine 
TZ e and IZ^ for 10 bins ranging from downward going 
(cos 9 — 1) to upward going (cos# = —1) neutrinos. 
We also allow for a simple, but more realistic, energy 
dependence for n(E), taken from additionally, we 
introduce some overlap of the bins. We compare our re- 
sults for the azimuthal dependence of the neutrinos to 
the dependence of the observed recoil electrons seen at 
Super-K, normalized to their no-oscillation Monte Carlo 
simulation, in Fig. |31 Though we do not model the re- 
coil electron, there is a strong correlation between the 
two processes for the high-energy events. The results are 
encouraging. For the lower energies, all the solutions pro- 
duce little electron neutrino oscillations as indicated by 
the data. However, there is a visible low-energy muon 
neutrino oscillation which is larger in the theory than in 
the data. An improved model of the atmospheric data 
is required to better understand this. Most importantly, 
the high-energy muon zenith-angle data is qualitatively 
similar to the results given by our model. 

In summary, within the model developed in Ref. [lfij ] 
of the neutrino oscillation data, we find that the region 
#13 < plays an important role in understanding the 



oscillation parameters for three-neutrino oscillations. As 
oscillation probabilities for negative and positive values 
of 6*i3 are not simply related, the analysis cannot be per- 
formed in terms of sin 2 #13. The work presented here is 
intended to motivate a more thorough and careful exam- 
ination of the 6*i3 < region of the parameter space. 
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